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Spin liquids occuring in 2D frustrated spin systems were initially assumed to appear at strongest 
frustration, but evidence grows that they more likely intervene at transitions between two different 
types of order. To identify if this is more general, we here analyze a generalization of the spatially 
anisotropic triangular lattice (SATL) with antiferromagnetic Heisenberg interactions, the spatially 
eompletely anisotropic triangular lattice (SCATL). Using Takahashi's modified spin- wave theory, 
complemented by exact diagonalizations, we find indications that indeed different kinds of order 
are always separated by disordered phases. Our results further suggest that two gapped non- 
magnetic phases, identified as distinct in the SATL, are actually continuously connected via the 
additional anisotropy of the SCATL. Finally, measurements on several materials found magnetic 
long-range order where calculations on the SATL predict disordered behavior. Our results suggest a 
simple explanation through the additional anisotropy of the SCATL, which locates the corresponding 
parameter values in ordered phases. The studied model might therefore not only yield fundamental 
insight into quantum disordered phases, but should also be relevant for experiments on the quest 
for spin liquids. 

PACS numbers: 75.10. Jm,75.10.Kt,75.30.Ds,75.30.Kz 



I. INTRODUCTION 

Understanding magnetically disordered quantum ma- 
terials is of fundamental interest, e.g., for layered mag- 
netic insulators/metals in which magnetism is disrupted 
by charge doping, leading to dramatic phenomena such 
as high-temperature superconductivity [1-3J. Also, dis- 
ordered quantum phases can have excitations which are 
fractionalized even in two dimensions [4J. However, clas- 
sical order is typically quite resilient, especially in two 
or three dimensions [5HS]- To disrupt the classical order 
and reach quantum-disordered phases like valence-bond 
solids or resonating valence-bond states, it is assumed 
that frustration could be a crucial ingredient [4l[9]. Here, 
we want to investigate where such quantum-disordered 
phases appear in a two-dimensional, antiferromagnetic 
(AFM) Heisenberg model with highly tunable frustra- 
tion. 

The Heisenberg Hamiltonian describes a large variety 
of magnetic materials. It reads 

where Si is a Heisenberg spin-S* operator at site i (here 
we are interested in the extreme quantum limit S = 1/2). 
We consider a triangular geometry where the nearest- 
neighbor (NN) couplings Jij along all three lattice direc- 
tions are different, the spatially completely anisotropic 
triangular lattice (SCATL), see left side of Fig. [l] For 
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Figure 1: Geometry of the SCATL. The spins (gray 
bullets) are coupled to NNs along the lattice vectors ti,2,3 by 
the couplings J-n = J, J1-2 = J\ and J-rg = J" ^ which can all 
be mutually different. The right part of the Figure shows the 
associated square lattice. The shown geometry is the one used 
in the ED of the 15-site system, chosen for maximal symmetry 
between all three couplings. 



simplicity, we will work throughout this paper in the asso- 
ciated square lattice (right side of Fig. [T]), where the vec- 
tors connecting NN sites are ti = (1, Ij? T2 = (0, 1), and 
Ta = (-1,0), and define J^^ = J, J^^ = J', and J^^ = J". 
This model generalizes the spatially anisotropic triangu- 
lar lattice (SATL), where two of the couplings are equal. 
Similar to the SATL, the SCATL can be tuned between 
vanishing and strong frustration. 

This specific model may be relevant to understand the 
spin liquids occurring in related systems. E.g., in a pre- 
vious work on the SATL and the square-lattice J1J2J3 
model (which has NN, next-NN, and next-next-NN in- 
teractions) [lOj, we found that spiral long-range order 
(LRO) is never connected directly to collinear order. In- 
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stead the system seems to pass through a quantum disor- 
dered region. In this article, we want to examine if such 
a behavior remains true in the more general SCATL, to 
further corroborate if it is a general characteristic of two- 
dimensional frustrated spin systems. 

A related advantage of the chosen geometry is the pos- 
sibility, given by the additional anisotropy, to approach 
the putative non-magnetic phases from different angles, 
possibly revealing crucial information not only about 
their location in parameter space, but also about their na- 
ture. Indeed, the results presented in this paper suggest 
that two gapped non-magentic regions, previously identi- 
fied as two distinct phases in the SATL, might actually be 
continuously connected via the additional anisotropy of 
the SCATL. Studying the persistence and characteristics 
of the putative quantum-disordered phases with respect 
to this additional anisotropy is the first main aim of the 
present paper. 

The second main aim is related to experimental find- 
ings in magnetic materials. While the SATL has found 
considerable attention in recent years, to our knowl- 
edge the ground-state phase diagram of the Heisen- 
berg SCATL has never been thoroughly investigated. 
Recent first-principles calculations, however, show that 
some magnetic materials, such as the organic salts 
Me4-nEtnPn[Pd(dmit)2]2 (abbreviated Pn - n) [27J, 
TMTTF [28j, or BaAg2Cu[V04]2 [29j, which are weh 
described by weakly-coupled 2D triangular lattices, can 
have considerable anisotropics between all three intra- 
plane couplings [66j. Typically, to locate the material 
within the well-studied SATL model, the two closer bond 
strengths are averaged. However, this places materials 
such as Sb-0 and As-2, which are experimentally found 
to be AFM ordered, into a region of the phase diagram, 
where according to many theoretical studies fiO\ l3Qti33] 
no LRO should exist. The second aim of this paper, 
therefore, is to show that the additional anisotropy be- 
tween the couplings could naturally explain this discrep- 
ancy, since it shifts the parameter values corresponding to 
these materials into an ordered phase. This also suggests 
that the non-magnetic state is quite sensitive to this ad- 
ditional anisotropy, which therefore has to be taken into 
account when interpreting experiments. 

In this work, we investigate the S = 1/2 Heisenberg 
AFM SCATL, Eq. ([T]), within Takahashi's modified spin- 
wave theory (MSWT) [34J, supplemented with the op- 
timization of the ordering vector [10^, ^35\. Previously 
[TOl [35] , we have shown that this improves significantly 
over conventional spin-wave theory (as well as over con- 
ventional MSWT), as it allows to account for the dra- 
matic quantum corrections to the type of order appear- 
ing in frustrated quantum antiferromagnets. Further, the 
breakdown of the theory provides a strong signal that the 
true ground state might be quantum disordered; hence, 
this method serves to efficiently find candidate models for 
spin-liquid behavior. While the main focus of this article 
is on the MSWT results, we complement them with ex- 
act diagonalization (ED) of small clusters. The SCATL 



with XY interactions, motivated by recent experiments 
with frustrated bosonic atoms in optical lattices, will be 
treated in a similar way in the following article |36| . 

The rest of this paper is organized as follows. First, 
to understand which effects can be expected in the quan- 
tum SCATL, we discuss the phase diagram of its classi- 
cal counterpart (Sec. I A) and briefly summarize known 
results from its well- studied limiting case, the SATL 



(Sec. IB). Sec. [TT] contains the main results of our pa- 
per, namely the discussion of the quantum-mechanical 
ground-state phase diagram of the SCATL, including var- 
ious observables from MSWT and ED as well as, for a 
possible comparison to experiment, the spin- wave disper- 
sions at selected points of the phase diagram. We dele- 
gate the technical details of the MSWT to the Appendix. 
Sec. |III[ finally, provides some conclusions. 



A. Classical phase diagram 

In this section, we discuss the classical phase diagram 
of the SCATL, which can serve as a guide to what ordered 
phases are to be expected, and which allows to appreciate 
the changes brought about by quantum fluctuations. 

To obtain the classical solution, we replace the Heisen- 
berg spins in Eq. ([T]) by classical rotors, which - without 
loss of generality - lie in the xy-pldiue [67\. The order- 
ing vector Q^^ = {Q^jQ*^) is the /c-vector which mini- 
mizes the Fourier transform of the coupling strengths. It 
fixes the direction of each spin (up to a global phase) as 
Si = S (cos(Q^^ • n), sin(Q^^ • r^)). We find 
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The classical phase diagram of the SCATL, plotted 
in Fig. |2j contains several Neel-ordered phases and an 
extended spiral-ordered phase. The Neel phases spread 
around the square-lattice limits {J' I J" I J) = (1,0) 
with = (0,7r), {J' I J, J" I J) = (0,1) with = (7r,0), 
and J' I J, J" I J » 1 with Q^^ = (7r,7r). The spiral 
phase, with continuously varying ordering vector, con- 
nects smoothly to the Neel phases, and occupies the ex- 
tended region between them. In particular, it extends 
all the way to J'/ J = J" j J = [and, symmetrically, to 
( J7J = 1, J'VJ ^ oo) and (J'VJ = 1, J7J ^ oo)], where 
the system decouples into an ensemble of ID chains. 



B. Known results in limiting cases 

In this section, we discuss well-known limiting cases 
of the quantum SCATL, including results on the SATL. 
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Figure 2: Classical phase diagram of the SCATL: The 
ordering vector evidences tiiree Neel-ordered phases, an ex- 
tended spiral-ordered phase, and limits where the system de- 
couples into an ensemble of independent chains, as indicated 
by the labels in the left panel. The thick red lines denote tran- 
sitions between different kinds of order, and along the dashed 
black lines the system is in SATL limit. 



This helps us to assess which phases and quantum effects 
are to be expected in the phase diagram of the SCATL. 

For J' I J, J" I J » 1, {J' I J, J" I J) = (1,0), and 
{J' I J" I J) = (0,1), one recovers the square lattice 
limit. Here, Neel order persists also in the quantum 
case [21]. Similarly, in the isotropic triangular lattice, 
J' = J'' = J, spiral LRO survives quantum fluctuations 
p7j. The limits (J' = J'' = 0), (J' ^ oo with J'' = const), 
and (J'' oo with J' = const) correspond to ensembles 
of decoupled, critical Heisenberg chains, with algebraic 
correlations along individual chains but no correlations 
between them. 

For J' = J'' = a J (or, equivalently, J' = J or J'' = J), 
one recovers the SATL, which is realized in a variety of 
S = 1/2 compounds, e.g., CS2CUCI4 [38j and /^-(BEDT- 
TTF)2Cu2(CN)3 [Ml SO]. The model may display spin- 
liquid phases, although their extent and nature is still 
under intensive debate [TOl HlHi6]. In the rest 

of this section, we review the main features of the SATL 
phase diagram as found in the literature, proceeding from 
large to small a = J' j J = J" j J (for comparison. Fig. [s] 
reproduces the MSWT phase diagram from Ref. [lOj). 

It is commonly accepted that order-by-disorder effects 
due to quantum fluctuations stabilize the Neel phase 
considerably over the classical model, moving the point 
where Neel order disappears downwards from the classi- 
cal value a = 2 to values between a ^ 1.1 and 1.67, de- 
pending on the method used [iQl [33l SB SSI SJ, '48J . Fur- 
ther, several methods predict that quantum fluctuations 
spread the transition point between the Neel and the spi- 
ral phase into a quantum-disordered phase [TO', "41] . 
In the following, we term this predicted disordered region 
"large-Qf quantum-disordered region" (large-a QDR). One 
of the main aims of this article is to study if it is a gen- 
eral feature of frustrated quantum antiferromagnets that 
a quantum-disordered phase intervenes in transitions be- 
tween commensurate and incommensurate order. 



Similar behavior has been found in a variety of quan- 
tum spin models, including Ji J2 Js-models on the square 
lattice [TDHH], and frustrated honeycomb models with 
Heisenberg [23^-25] or XY interactions [26j. In fact, since 
quantum phase transitions are driven by quantum fluctu- 
ations, one might expect that - if anywhere - a complete 
restructuring of the ground state in favour of a quantum 
mechanical configuration may occur preferably close to 
a quantum critical point. It, hence, seems plausible that 
at such points quantum fluctuations are most effective in 
disrupting classical order. Indeed, a similar effect occurs 
in classical statistical physics. Assume that there is a 
transition between a commensurate and an incommen- 
surate phase which both show LRO. As a first notable 
thing, close to this transition the thermal phase transi- 
tion to a disordered state will typically happen at lower 
temperature than far away from it. Beyond the thermal 
phase transition, the disordered phase will show short- 
range order of the type corresponding to the adjacent 
long-range ordered phase. The transition between the 
two different kinds of short-range order is called a dis- 
order point [1SH52] . Interestingly, the correlation length 
associated with the two kinds of short-range order has a 
minimum just at this point. Hence - similar to what is 
found in the quantum models - thermal fluctuations tend 
to suppress order most effectively at a commensurate- 
incommensurate transition. 

At the low-a side of the spiral phase, previous works 
predict a disordered phase, which could appear for as 
large a as ^ 0.8 - 0.9 [30-33J. Our previous MSWT re- 
sults suggest a ^ 0.65 [lOj. In the following, we term 
this predicted disordered region "small-a QDP." It may 
be associated to a spread of the gapless spin liquid of the 
isolated chains {J' = 0) to finite coupling [3TH33j , possibly 
followed by a gapped spin liquid [10\ [5T| [3^ . This double 
nature of the disordered region is still under debate, since 
some works only find a gapless spin liquid [33j. Consis- 
tent between these methods is the prediction that quan- 
tum fluctuations disrupt ordering tendencies between the 
chains even for relatively large inter-chain couplings. But 
consent about the physics in this region seems far from 
reached. For example, recent DMRG studies entirely 
question the existence of the small-a spin liquid(s) [45]. 
And a recent renormalization-group analysis [46j found 
collinear AFM long-range order in the region a < 0.3 (see 
also [44J), and, above that value, spiral order. The weak- 
ness of the spiral order leaves, however, the possibility 
that the true quantum ground state hosts a disordered 
phase in the parameter range 0.3-0.5. 



II. QUANTUM-MECHANICAL PHASE 
DIAGRAM 

From the discussion of the classical phase diagram and 
the limiting cases, we have the necessary background to 
tackle the quantum-mechanical ground-state phase dia- 
gram of the Heisenberg SCATL. To compute it, we use 
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Figure 3: The MSWT quantum phase diagram of the 

SATL (from Ref. [10 ) contains Neel order (which is consid- 
erably more stable than in the classical model), spiral order 
(which is destabilized by quantum fluctuations), and two pu- 
tative spin-liquid (SL) phases. These are found through the 
breakdown of the theory and a disappearing spin stiffness, 
which indicates a gapped disordered phase. In the purely ID 
limit, MSWT recovers convergence and produces the ID crit- 
ical state. We include sketches of classical states, where blue 
arrows indicate the directions of the classical rotors, namely, 
the ID state at J' /J = 0, the spiral state at J' /J = 1, and the 
2D-Neel state at a > 2, . 



the MSWT supplemented with ordering- vector optimiza- 
tion, working directly in the thermodynamic limit. Since 
this method is described in detail in our previous ar- 
ticles [10,. ^5j, we delegate the technical aspects to the 
Appendix, and only summarize here the main idea. The 
starting point is a classical state, which one dresses with 
quantum fluctuations in a second-order spin- wave expan- 
sion. This yields a bosonic Hamiltonian, the ground 
state of which is found self-consistently by minimizing its 
mean- field free energy. For this, quartic terms, i.e., in- 
teractions between spin waves, are decoupled via Wick's 
theorem. Additionally, we employ Takahashi's modifica- 
tion of vanishing magnetization. This constricts the aver- 
age number of spin- wave excitations to a physical value, 
in contrast to conventional spin-wave theory, where the 
spin- wave excitations grow completely unchecked. This 
modification has proven a crucial improvement to de- 
scribe low- dimensional systems with weak order tenden- 
cies. 

Typically, one uses the classical ground state as the 
reference state. However, in many models quantum fluc- 
tuations considerably shift the type of predominant or- 
der. Therefore, we find the ordering vector giving the 
best classical reference state by including it in the self- 
consistent optimization. This has proven crucial to cap- 
ture, e.g., the stabilization of the Neel phase by quantum 
fluctuations. As has been proposed in Refs. [I0l|35j, the 
breakdown of the theory strongly suggests that at mean- 
field level no semi-classical reference state yields a good 
description of the quantum ground state. This is then 
interpreted as an indication of non-magnetic behavior in 
the true ground state. This will be an important aspect 
for the interpretation of the quantum phase diagram. 

We compare these MSWT results to exact diagonaliza- 
tion (ED) of a 15-site lattice, as depicted in Fig. [l] The 



geometry is chosen for its symmetry between J, J', J'' 
bonds. It is important to leave the boundaries open to 
allow for incommensurate ordering vectors. 



A. MSWT and ED results — ordering vector and 
order parameter 

In this section, we give a first overview over the phase 
diagram, as obtained from the ordering vector Q and the 
order parameter M, followed in the next two sections by 
more detailed analyses. In MSWT, ordering vector and 
order parameter are direct results of the optimization 
[see Appendix, Eqs. (A8) and (All)]. In ED, they can 
be extracted from the static structure factor. 
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Its peak lies at the ordering vector Q^^, and the square 
root of its height, = M^^, approaches in the 

thermodynamic limit the order parameter M. 

As seen in the MSWT and ED ordering vectors, pre- 
sented in Fig. [4j quantum fluctuations stabilize the Neel 
phases compared to the classical case, as already observed 
in the SATL (Sec. IB). In the central region around 
J' J'' f« J, a broad range of incommensurate order- 
ing vectors indicates spiral order. The finite MSWT or- 
der parameter (Fig. [sj left panel) shows that in these 
phases indeed LRO survives quantum fluctuations. (Note 
that the self-consistent MSWT calculations become rela- 
tively unstable for small order parameters, which results 
in ragged phase boundaries.) 

In the Neel phases, the ED order parameter (Fig. [sj 
right panel) is maximal, giving support to the assumption 
that here LRO persists. However, it is much smaller in 
the spiral phase than the MSWT value, a discrepancy 
already found in the SATL [lOj. (This could be due to 
third-order corrections to the spin- wave expansion which 
our approach neglects, and which can become important 
in spiral configurations [53j.) 

Between the ordered phases, we find a broad region 
where MSWT breaks down, indicating as usual fiOi [35] 
that these regions do not allow a description in terms 
of an ordered, semi-classical state. Therefore, it ap- 
pears that it is a quite universal feature of frustrated 
quantum antiferromagnets that spiral- and collinearly- 
ordered phases are always separated by quantum disor- 
dered phases. This is the first main result of this paper. 

The strong decrease of the MSWT and ED order pa- 
rameters (Fig.|5| upon approaching this region gives sup- 
port to this interpretation (which we will further corrob- 
orate in the next two sections). Note also that both the 
ED and MSWT order parameter seem to disappear more 
smoothly when approaching the putative ID-like QDR 
(consider, e.g., in the range 2 < J' /J < 3, J''/ J ^ 1~). 
Upon approaching the putative large-a QDR dividing 
spiral from Neel LRO, on the other hand, for ED, the or- 
der parameter decreases sharply (consider, e.g., the line 
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Figure 4: Quantum- mechanical phase diagram of the 
SCATL, ordering vector. Upper row: MSWT data. 

Quantum fluctuations stabilize the Neel phase. Around J' ^ 
J" J, a part of the classical spiral phase survives quan- 
tum fluctuations (labels in the upper left panel). The solid 
symbols (upper right panel) denote some experimental mate- 
rials (see Table (For clarity, we show only symbols in the 
lower right part of the figure, excluding points symmetric un- 
der exchange of J, J' ^ and J" .) A: magnetically disordered, 
• : charge ordered, ■: AFM LRO. Note especially the two ■ 
at {f I J, J" I J) = (1.44,0.84) and (f/J,J''/J) = (1.36,0.82) 
marking the materials As-2 and Sb-0, which lie well inside a 
Neel ordered phase. Neglecting the asymmetry between the 
couplings would put As-2 into the supposedly disordered re- 
gion and Sb-0 just at its boundary (□). Lower row: ED 
data for N = 15 sites. Already for this small system, it can 
be appreciated that (compared to the classical case) the Neel 
phase grows at the expense of spiral order. 




Figure 5: Quantum- mechanical phase diagram, or- 
der parameter. ED results qualitatively confirm MSWT. 
In particular, the order parameter for both methods decreases 
rapidly upon approaching the MSWT breakdown regions. 



material 


/ -rf 1 T -rff 1 T\ 
(J /J, J /J) 


state 


N-3 


(10, 9.1) 


AFM 


P-2 


(2.63, 1.89) 


AFM 


Sb-2 


(2.08, 1.79) 


CO 


Sb-1 


(1.72, 1.49) 


SL 


K-CN 


(1.41, 1.41) 


SL 


P-1 


(1.32, 1.32) 


VBS 


Cs 


(1.16, 1.16) 


CO 


Sb-0 


(0.74, 0.60) 


AFM 


Cs2CuBr4 


(0.74, 0.74) 


AFM 


As-2 


(0.69, 0.58) 


AFM 


CssCuCU 


(0.34, 0.34) 


RVB 



Table I: Some relevant materials for which the ground 
state has been measured in experiment, together with the cou- 
pling strengths, and the state they are found to be in (from 
Refs. [271 [56l |57] and references therein). AFM: antiferro- 
magnetic LRO, CO: charge ordered, SL: spin liquid, VBS: 
valence-bond solid, RVB: resonating valence-bond state. 



J7J = 1, J'VJ ^1"). Here, for MSWT, the break- 
down occurs abruptly at finite order parameters. This 
could point at a difference in the type of phase transi- 
tion upon approaching the large-a QDR and the non- 
magnetic phase at the decoupled-chains limit. 

The second main result of our paper concerns ex- 
perimental measurements of ground-state behavior of 
some materials, taken from Ref. [27\ (see also the re- 
views [ST, 'SS]), as well as from Refs. [56j (CS2CUCI4) 
and [57J (Cs2CuBr4). For reference, they are presented in 
Table |Tj and included as solid symbols in the upper right 
panel of Fig. [4j We mark magnetically disordered ma- 
terials (spin liquids, resonating valence-bond states, or 
valence-bond solids) with triangles, charge-ordered ma- 
terials with bullets, and AFM-ordered materials with 
squares. Including the full anisotropy of the triangular 
lattice, all AFM ordered materials lie within the ordered 
phases from MSWT [68j. In particular, the AFM ordered 
materials As-2 and Sb-0 [at (J'/J.^/J) = (1.44,0.84) 
and {J' I J, J" I J) = (1.36,0.82)] lie well inside a Neel 
ordered phase. If one neglects the anisotropy between 
J' and J'', taking as usual the mean of both couplings, 
they would lie at the position of the empty squares at 
{J' I J, J" I J) = (1.57,1) [equivalent to {J'lJ,J"lJ) = 
(0.64,0.64)] and {J' /J, J'' /J) = (1.50,1) [equivalent to 
{J' jJ^ J" I J) = (0.67,0.67)] - inside a phase where many 
methods [3QH33] predict disorder; specifically, within 
MSWT, symmetrizing the couplings puts Sb-0 just at 
the border to the breakdown region (which should be 
a lower limit for a disordered phase in the true ground 
state) and As-2 within it. The appearance of AFM Neel 
LRO in these experiments may find, therefore, a simple 
explanation in the full anisotropy of the SCATL. This 
second main result of our paper shows how crucial the 
full anisotropy is for the interpretation of experimental 
data. The rest of this article is devoted to fleshing these 
main findings out. 
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Figure 6: Upper panel: The partial spin stiffness 

decreases upon approaching the MSWT breakdown region, 
suggesting the disruption of magnetic LRO. Lower panels: 
The spin- wave velocities perpendicular to the dominating 
coupling strength soften in the ID limits. Differences in the 
spin-wave velocities might allow to measure the anisotropy 
of the SCATL. All quantities are normalized to the coupling 
strengths 1 -\- J' -\- J" . 



B. Supporting observables from MSWT 
stiffness and spin-wave velocities 



spin 



In Refs. [T0l[35], the spin-stiffness tensor, which char- 
acterizes the stiffness of the magnetic order under change 
of the ordering vector, has proven a valuable consistency 
check of our MSWT calculations. Its components are 



dQadQfs ' 



(4) 



where is the free energy. Even if the order parameter 
is finite, a small spin stiffness suggests that further quan- 
tum fluctuations than taken into account within MSWT 
could disrupt the remaining order [69 j. 

Since for our purposes an upper bound for the spin 
stiffness is sufficient, we take the partial derivative in 
Eq. Q. The exact spin stiffness can be computed via 
the total derivative. To this, within the self-consistent 
MSWT calculations, one first has to find the optimal 
ordering vector. Then, one reruns the self-consistent 
MSWT equations for several fixed, slightly non-optimal 
ordering vectors, yielding slightly larger energies. The 
spin stiffness can be derived by fitting a quadratic form 
to the resulting energy landscape. In the self-consistent 



iteration, the mean fields characterizing the MSWT state 
can adjust to a changed ordering vector. This effect is not 
taken into account in the partial derivative, which hence 
provides an upper bound to the total spin stiffness. We 
find that it suffices to extract the location of disordered 
phases, but it may yield wrong results about their nature. 
In particular, we found in the SATL [lOj that, upon ap- 
proaching the putative small-a QDR, not only the total 
inter-chain, but also the total intra-chain spin stiffness 
decreases strongly. Since such a behavior is not consis- 
tent with algebraic correlations along the chains, this can 
be interpreted as an indication of a gapped quantum- 
disordered state. The partial spin stiffness computed in 
Ref. [lOj, on the other hand, only vanishes in the inter- 
chain direction. Hence, it may not be able to distinguish 
gapped from gapless spin liquids. However, it still seems 
to adequately capture the location of disordered regions. 

In Fig.[6l upper panel, we show the determinant of the 
spin-stiffness tensor, det(p), normalized to the coupling 
strengths 1 + J' + As we should expect [58j, det(p) 
decreases upon approaching the phase transitions, espe- 
cially from the Neel-ordered side. At large J' (J''), this 
decrease is due to a softening of the stiffness in x (y) 
direction, and at small {J' jJ^ J" I J) the direction per- 
pendicular to Ti (as has also been found in Ref. [lOj). 

Another indicator for approaching disordered phases 
is given by the spin- wave velocities Vx^y^ which can be 
connected to the spin stiffness via the susceptibility [51] . 
Since the spin-wave velocities are defined as the leading 
order of an expansion of the spin-wave dispersion rela- 
tion, Eq. (A7), around small i.e.. 



lim Ukjkx 



lim oJk/ky 



ky=0 



(5a) 



(5b) 



they can be measured directly from the spin- wave disper- 
sion, allowing an experimental check of our findings. 

As seen in Fig. |6j lower panels, close to the ID break- 
down region, they, too, soften in the direction perpen- 
dicular to the dominating coupling. On the other hand, 
when approaching the putative large-a QDR dividing the 
spiral from the Neel phase, both spin-wave velocities re- 
main finite. This is another (besides the different be- 
havior of the order parameter) indication that the large- 
a QDR could be qualitatively different from the non- 
magnetic phase found in the limit of decoupled chains. 



Supporting observables from ED — energy 
derivative, gap, and chiral correlations 



The ED observables investigated in Sec. II A allowed to 
interpret the predominant ordering behavior, but did not 
yield clear evidence if within ED really quantum phase 
transitions exist, and if yes, where. The second deriva- 
tives of the ED ground-state energy per spin, plotted in 



Lid ^E^^/a J'^) Lid ^E^^/d J"^) 



Figure 7: Second derivative of ED ground-state en- 
ergy per spin for N = 15. For clarity, we plot the logarithm 
after a shift to values larger one, L(d^E^^/dJ^^), where 
L{x) = log(l + max(x) - x), and = J' or J". Strong peaks 
mark the transitions from the Neel phases. An additional 
peak around {J'^J") = (1,1) might be an indication of an 
additional phase, separating the Neel phases from the spiral 
one. 

Fig. [7| can provide such an indicator. In the thermody- 
namic limit, it diverges at a quantum phase transition. 

Indeed, there are clear peaks at lines similar to where in 
MSWT the Neel order breaks down. Also, a peak appears 
around {J' ^J") = (1,1). This might be a precursor of a 
quantum phase transition away from the spiral state, and 
to an intermediate phase possibly to the QDR which is 
supposed to exist in this system. 

We get further support for this phase diagram from the 
ED energy gap between ground and first excited state. 
Fig. [Sj In the well-known limiting cases of the SCATL, 
it behaves as expected: There is no singlet gap close to 
the decoupled-chains limits, since the system is then in 
a critical phase. In the Neel ordered phases, there is a 
large gap which separates the ground state from closely- 
spaced excitations, which in larger lattices become the 
spin waves, collapsing slowly towards the ground state 
[60j. This is consistent with the considerable size depen- 
dence found for these parameter regions, as can be seen 
in the right panels of Fig. [8) where we plot cuts of the 
gap per spin, A£^™, at fixed J'^/J = 1, 2, 3 for triangular 
systems similar to the one in Fig. [TJwith N = 6^ 10, 15. 

On the contrary, there is no gap in the spiral phase, 
because there are two degenerate ground states with op- 
posite chirality [70j. We find that this vanishing of the 
gap depends strongly on the system geometry, but it oc- 
curs consistently for all triangular systems considered. 

Interestingly, the gapless spiral phase is surrounded by 
a region where the gap attains considerable values. The 
very small dependence on system size for this parame- 
ter region indicates that this is stable towards the ther- 
modynamic limit. A finite gap is not consistent with a 
spiral-ordered phase. On the other hand, the predom- 
inant order in this region is at incommensurate wave- 
vectors. Hence, the finite gap is clearly not due to square- 
lattice Neel physics. Optimistically, these findings could 
therefore be interpreted as the precursors of a gapped 
QDR. This gapped region completely encircles the spi- 




Figure 8: Left panel: The ED gap per spin lends support 
to the MSWT phase diagram: A finite gap separates in the 
Neel phases spin- wave excitations from the ground state. In 
the spiral phase, the ground state is doubly degenerate due to 
the ambiguity in choice of chirality. The finite gap surround- 
ing the degenerate region could be a precursor of a gapped, 
disordered phase. At the quantum phase transitions to the 
Neel phases, the gap closes again. Right panels: cuts at 
fixed J" /J - 1,2,3 for triangular systems with increasing N 
(from light to dark and thick to thin: 6,10,15). There is lit- 
tle size dependence in the central gapped phase {J" I J = 2,3 
with J' ^ J, as well as J^^ = J and J' / J > 1.5). Also, for 
J'^ = 3 J the transition points do not show any appreciable 
size dependence, while for J'' = 2 J the one around J' /J = 1.5 
does. 

ral phase, suggesting that the low- and large-a gapped 
QDRs found in the SATL could actually be continuously 
connected via the additional anisotropy of the SCATL. 

Upon approaching the Neel phases, the gap closes, in- 
dicating a quantum phase transition. 

Going back to Fig. [5j when comparing MSWT and ED 
the lateral extent of the putative non-magnetic phases 
is different. Scanning along J' = J, within MSWT it 
is smallest around J'' = 2J, while for ED it decreases 
monotonously with increasing J" . As the right panels in 
Fig. |8] indicate, this discrepancy could be due to finite- 
size effects. Indeed, we find in ED that for J'' /J = 2 the 
transition point at J^/J = 1.5 shows an appreciable size 
dependence, while the transition points for J^'/J = 3 do 
not. Therefore, the lateral extent of the putative QDRs 
at around J^/J = 2 could decrease with making the 
MSWT and ED pictures consistent. 

From the gap, it seems that there is support for an 
extended gapped phase separating spiral and Neel LRO. 
Still, it would be desirable to exclude for this region spiral 
LRO in the thermodynamic limit. For this, we now study 
where chiral correlations persist. The vector chirality is 
defined as 

2 

Kij^k = {Si X Sj + Sj X Sk + Skx Si)^ , (6) 

where the sites /c} are located counter-clockwise on 
a triangle. For the small systems used in our ED, we 
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N=12 



7V=18 



A^=7 



Figure 9: Chiral correlations from ED. Upper left: 

Already for small systems (N = 15), chiral correlations are 
appreciably smaller in the Neel phases than in the rest of 
the phase diagram. Upper right: Extrapolating to large 
lattices, chiral LRO apparently only survives in a small region 
around {J' /J, J" I J) = (1? I), lending support to an extended 
non-magnetic phase surrounding a spiral phase. Below: The 
geometries used in the extrapolation are chosen for symmetry 
upon rotation by 60° and equal number of J, J\ and J" 
bonds. 



generalize the chiral correlations [61 J to 
4 



(7) 



Here, the sum a runs over all triangles, while c runs only 
over the central ones to reduce boundary effects. The 
factors 5a, c weight Ka^c with a + (-) sign if the triangle 
points upwards (downwards). The prefactor, where A^a 
is the number of summands, is chosen such that the chiral 
correlation has the same theoretical maximum of | as the 
usual definition for large lattices |61| . 

As can be seen from the ED results of the A/" = 15 
lattice (Fig. |9j left panel), the chiral correlations are 
relatively small in the Neel phases and largest in the 
spiral phase around {J' j J" j J) ^ (1,1). However, at 
this lattice size, there are still appreciable chiral cor- 
relations in the rest of the parameter regime. In par- 
ticular, in the ID limit, they are only a little smaller 
than in the spiral phase. Therefore, we also plot in 
Fig. [9j right panel, an extrapolation to large lattices by 
^f_{N) = ^_(7V = oo) + + 1^ + where we use the 

known form for the leading finite-size behavior [62j but 
also include subleading corrections due to the small sys- 
tems under consideration (our data comes from lattices 
with N = 7^ 10, 12, 15, 18, all chosen to have the same 
number of J, J\ and J'' bonds, as sketched at the bot- 
tom of Fig.^. This shows a clear trend, namely that the 



chiral correlations only survive in a small region around 
{J' I J" I J) = (1,1), roughly where the vanishing gap 
indicated the spiral phase [71j. This would mean that 
outside this region there is no spiral LRO. 

With this, we have several independent observations 
from ED indicating the existence of a magnetically dis- 
ordered phase surrounding the spiral phase: the in- 
crease of the gap when leaving the central region around 
{J' I J" I J) = (1,1) and the disappearance of chiral 
LRO for large lattices both suggest that there is no spiral 
LRO in this region. On the other hand, the predominant 
order is at incommensurate ordering vectors, indicating 
that this phase is also not Neel ordered. Therefore, it 
seems natural to assume that this region could host a 
non-magnetic phase, possibly gapped far away from the 
ID limit and gapless close to it, consistent with MSWT. 



D. MSWT spin-wave dispersion relations 

Finally, to connect to experiment, we provide the spin- 
wave dispersion relations cj^, as computed from MSWT, 
Eq. (A7). In Fig. 10, we show parameters corresponding 



to a point from the spiral phase and the magnetically or- 
dered materials listed in Table [ij We also provide (where 
applicable) a comparison to the dispersion relation which 
would result if two of the couplings were equal. These 
comparisons can be seen more quantitatively in the cuts 
(c.i-iii) shown in the lowest row of Fig. 10 For the point 
from the spiral phase (a.i), the symmetrization (b.i) does 
not significantly change the dispersion relation, but for 
P-2 and, especially, for Sb-0, the differences are consider- 
able. The latter in particular changes even qualitatively 
since a symmetrization would put it instead of into a 
Neel phase into a spiral phase. These differences seem 
significant enough to be measurable in experiment. Such 
a measurement could allow to quantify the actual mag- 
nitude of coupling anisotropics. 



III. CONCLUSION 

In conclusion, we have provided a thorough analysis of 
the ground-state phase diagram of the quantum Heisen- 
berg SCATL. Using various observables from modified 
spin-wave theory supplemented with ordering- vector op- 
timization, and supported by exact diagonalization data, 
we have found that quantum fluctuations stabilize Neel 
order with respect to the classical phase diagram. Fur- 
ther, they reduce the extent of the spiral phase, which 
seems to be entirely surrounded by a quantum disordered 
region. This result, which constitutes our first main find- 
ing, is supported by the breakdown of MSWT, together 
with the strong decrease of the order parameter and the 
spin stiffness. While MSWT cannot be applied to study- 
ing this region, the fact that no semiclassical reference 
state describable by an ordering vector yields a stable 
solution is highly suggestive of a magnetically-disordered 
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Figure 10: Spin- wave dispersion from MSWT, normalized to 1 + J' + J" . Top row: parameters corresponding to 
a completely anisotropic point in the spiral phase (a.i), and the magnetically ordered materials cited in Ref. [27 (a.ii-iv). 
Middle row: symmetrizing the two closer couplings can change the dispersion relations (b.i-iii). (b.iv): dispersion relation 
for Cs2CuBr4. Bottom row: Cuts along the path indicated in (b.i) allow to more quantitatively compare the completely 
anisotropic dispersion relation taken from row (a) (thick line) to the corresponding symmetrized one from row (b) (thin line). 



phase of purely quantum origin. Hence, our results out- 
line a very promising candidate region for such exotic 
states, meriting further research with more sophisticated 
theoretical methods or experimental setups. 

The possible existence of quantum-disordered phases 
is further corroborated by ED data, where a finite gap 
and a vanishing chiral correlation make spiral LRO seem 
unlikely, while the location of the structure-factor peak 
at incommensurate wave vectors seems to preclude Neel 
LRO. Also, the strong decrease of the ED structure- factor 
peak appears to support this interpretation. Further, we 
found some indications that very close to the ID limit the 
transition from the Neel phase to the putative disordered 
region could be qualitatively different from what happens 
at larger a. 

A complete encircling of the spiral phase by disor- 
dered phases could naturally explain the succession of 
a gapped and a gapless non-magnetic phase at the low- 
a limit of the SATL. The gapless quantum-disordered 
phase would be continuously connected to the limit 
of decoupled chains, while the additional anisotropy 
of the SCATL would adiabatically connect the gapped 
quantum-disordered phases at small and large a. There- 
fore, the additional anisotropy has great potential to de- 



liver new insights into the nature of these kind of phases. 

Our second main finding is connected to experimental 
results: measurements find magnetic LRO in materials 
which theoretical analyses on the SATL predict to be 
magnetically disordered. We show that this discrepancy 
finds a simple explanation in the additional anisotropy 
of the SCATL, which is neglected in the SATL. Taking 
it into account, we predict these material to lie in mag- 
netically ordered phases, in accordance to experiment. 
These findings show the importance of the complete lat- 
tice anisotropy for the explanation of recent experiments. 

Finally, we provided spin- wave dispersion relations, a 
comparison to which might allow to probe the additional 
anisotropy experimentally. 
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Appendix A: MSWT formalism 

In this Appendix, we shortly review the MSWT for 
Heisenberg antiferromagnets (see [lOj; for a fuh descrip- 
tion of the approach - as apphed to XY models - see [35j). 

A fundamental assumption of spin- wave theory is that 
the ground state has LRO with ordering vector Q. Hence, 
it is convenient to rotate the local reference system as 



5f 



- sin {Q-n) + cos {Q-n) , (Ala) 
cos{Q-ri)S^ + sm{Q-ri)S^ , (Alb) 
-Sf. (Ale) 



lei to the classical spin Si = (cos (Q • Vi) ^ sin (Q • r^) , 0). 
This defines the classical reference state. We do not make 
any assumption on the ordering vector Q. In particular, 
it may well differ from the one of the classical limit {Q^^). 

Spin waves around this reference state can be described 
by the Dyson-Maleev (DM) transformation [51], 
which maps the physical spins to interacting bosons. 



Si 



^(2S-alai)ai 



2Sal, 
S + aja,- 



(A2a) 

(A2b) 
(A2c) 



Then 6'/, which will be the quantization axis, lies paral- where = ± iS^. 



This transformation maps the Hamiltonian, Eq. Q, to the non-linear bosonic Hamiltonian 
H = ^ Jij{ [25' {alaj + a^a]) - ala^jajaj - ajaittia]] (1 + cos (Q • Tij)) 



+ [2S' (aja] + aittj) - aiO^-ajaj - a\aiaiaj\ (1 - cos (Q • Vij)) 
+ A[S'^ - S{a\ai + a\aj) + a\aia\aj]cos{Q-rij) + o{^ 



(A3) 



where {a\) destroys (creates) a DM boson at site i, 
and S is the length of the spin. Here, we neglected 
the kinematic constraint which restricts the DM-boson 
density n to the physical subspace n < 2S. Moreover, 
we dropped terms with six boson operators, which are 
of order 0[n/{2S)^] and are negligible for n/{2S) < 1. 
Using Wick's theorem [65j, and defining the correlators 
{alaj) = F (vij)- \5ij and {aiaj) = (a^^a]) = G {vij), the 
expectation value E = {H) can be written as 



E 



(A4) 



(1 + cos (Q • Tij)) 

(1 - COS (Q ■ Vij)) \ . 



5+--F(0)+F(r,,) 
5+i-F(0) + G(r,,) 



After Fourier transforming, = Ei cii e~^^''^^ , 

and a subsequent Bogoliubov transformation, 

cosh Ok ak - sinh^fctt^^, and = -sinh^^afc + 

cosh^fc a^^, we minimize the free energy under the con- 
straint of vanishing magnetization at each site, (a-a^) = 
5, which is known as Takahashi's modification [34J. This 
yields a set of self-consistent equations. 



tanh2^fc 



(A5) 



with 



^fc = ^ E -^ii (1 - cos (Q • Vij)) Gij e"'-'-'^ , (A6a) 



1 



(A6b) 



- (1 + cos {Q ■ nj)) (1 - e**'-"- )] - ^ , 

where /i is the Lagrange multiplier for Takahashi's con- 
straint. The spin- wave spectrum reads 



(A7) 



At T = 0, /i vanishes, which implies the disappearance of 
the gap ai k = that may exist for finite temperature. 
This is a necessary condition for magnetic LRO, and en- 
ables Bose condensation in the k = mode. Separating 
out its contribution, 

(4=o«'^=o) /N = {ak=oak=o) /N = M, (A8) 

which corresponds to the magnetic order parameter, one 
arrives at the zero-temperature equations 

F,j = M+^y^ COS {k • r,j) , (A9a) 

ktO 

= M + — ^ cos (fc • r,, ) , ( A9b) 
fc*0 ^fc 
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and the constraint of vanishing magnetization at each 
site becomes 



2 27V ujk 



(AlO) 



It is not a priori clear that the classical ordering vector 
Qc\ QQYYectly describes the LRO in the quantum system. 
To account for a competition between LRO at different 
ordering vectors Q, we extend the MSWT procedure by 
optimizing the free energy with respect to the ordering 
vector Q. This yields two additional equations which 
must be added to the set of self-consistent equations, 



dQy 



^--\T.J^3 sin (Q • n,)rl [F^^ + G^,] = . 



(Allb) 



d 



The values of Fij and Gjj can now be calculated by 



^ ,^^v>vv^^ ^"v^ I^J . . ^^^^^^^^^^ 

T^—Y. Jij sin (Q • rij)r'^j [Ff^ + G^J = , solving self-consistently Eqs. ( |A6HaTT] ). Through Wick's 



theorem the knowledge of the quantities Fij and Gij de- 
(Alla) termines the expectation value of any observable. 
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